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We show that the phase information contained in coherent phonon oscillations generated by a 
laser pulse in a multi- orbital superconductor can be used as an experimental tool to probe super¬ 
conductivity and pairing symmetries. The phase difference between the normal and superconducting 
states is proportional to the superconducting order parameter just below the superconducting tran¬ 
sition temperature, Tc. It also exhibits different behaviors for superconducting states with different 
pairing symmetries. In particular, if there is an orbital-dependent internal sign change state, the 
phase difference can have a discontinous jump below Tc. 


Introduction: The superconductivity in a multi-band 
electronic system can be extremly rich and complex. 
Many recently discovered correlated electron systems be¬ 
long to this category of multi-band superconductors. For 
example, iron-based superconductors discovered six years 
ago[T] have multiple Fermi surfaces and their bands near 
the Fermi level are attributed to all five d—orbitals. 
These materials exhibit a variety of intriguing proper¬ 
ties associated with all of the degrees of freedom includ¬ 
ing charge, orbital, spin and lattice [2], which can, in 
principle, lead to many possible novel superconducting 
states [3]. 

While theoretically, a multi band structure is a fertile 
ground for new physics, in experiments, it is still very 
difficult to detect them and determine their mechanisms 
because of the involvement of the multi-degrees of free¬ 
dom. Many experimental observations can be subject 
to multiple interpretations; for example, in iron-based 
superconductors [4], the interplay among electronic ne- 
maticity, magnetism and orbital ordering is still a subject 
of active research[5l |6]. The pairing symmetry of the su¬ 
perconducting state, arguably the most important prop¬ 
erty, is still controversial and highly debated [3]. While 
the magnitude of the superconducting order parameter 
can be directly probed by many experimental techniques, 
such as angle-resolved photoemission spectra (ARPES) 
and scanning tunneling microscopy(STM), there are few 
good direct probes to the phase distribution of the su¬ 
perconducting order parameter across their multiorbital 
Fermi surface. In particular, when the phase distribution 
is not enforced by the symmetry of the superconducting 
state, as the case stands in many theoretically proposed 
states in iron-based superconductors, the phase sensitive 
junction techniques [7] that determined the d—wave pair¬ 
ing symmetry in cuprates is also not applicable. 

Since the last couple of decades, ultrafast pump-probe 
spectroscopy has played an increasing role in probing the 
superconducting ground state, with the high Tc Cuprates 
having grabbed much of the attention [SHIT], along with 


a few experiments performed on multiorbital iron super¬ 
conductors [281136] as well. The primary focus of most 
of these experiments has been the measurement of re¬ 
laxation times that can be extracted from the behavior 
of the change in reflectivity ^ of the probe pulse as a 
function of the delay time 6 between the pump and the 
probe. From this, one can indirectly obtain information 
about the strength of the electron-phonon couplings and 
their anisotropies [TTJ |T2] , density of photoexcited quasi¬ 
particles mus], pseudo and superconducting gaps [25] , 
and signatures of the origin of the superconducting in¬ 
teraction However, even though coherent phonon 

oscillations in ultrafast experiments have been generated 
[Il[36| and studied [32l|34] for a while now, only a few ex¬ 
perimental works address the role of the superconducting 
phase on these oscillations and no theoretical background 
has been laid. 

In this Letter, we show that the phase of these coher¬ 
ent phonon oscillations contains useful information about 
the superconducting phase and its pairing symmetry; in 
particular, we show that the difference in the phase of 
the oscillations between the normal and superconducting 
state is proportional to the superconducting gap, and in 
certain scenarios, can help distinguish the sign change of 
superconducting orders on different bands. Thus, the co¬ 
herent phonons can act as a new experimental probe of 
superconducting symmetries. 

The coherent phonon amplitude mode with wave vec¬ 
tor q is described by the driven harmonic oscillator [37l [38] 

where Qq is the amplitude of the phonon mode, O is the 
frequency of the oscillator, (3 is the damping parameter 
and F{t) is the driving force. The solution to the above 
equation is given by 

Qq{t) = Ae~^^cos{^t + Fp/i) (2) 

where Q = and A is the amplitude of the 

oscillation which is proportional to the magnitude of the 
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driving force F. For simplicity, we will ignore any effect 
of damping. In such a case, the phase of the phonon 
oscillation Tph is given by j38] 


Tan{Tph) = 


Im{iF{-n)) 
Re{iF{-n)) ‘ 


(3) 


The driving force F{t) can be derived microscopically 
under reasonable approximations. Consider a general 
Hamiltonian that describes the physical processes in an 
ultrafast pump-probe experiment given by 


H = HeFHpF He-p + (4) 

where He, Hp, Hep and Hei{t) are electronic, phononic, 
electron-phonon coupling and electron-pulse interaction 
parts respectively [37]. In a superconducting state, the 
electronic part. He, is given by the general BCS form 

^/ccra/3c|,^Q,C/ccr/3 + ( 5 ) 

kaal3 ka 

where a and a are the orbital and spin index. We 
take the standard form for Hp = ^ 

Hep = J2kqaa' ^aa'Qqcl^Ck+ga' + h.c where Qg and Pg 
are the canonical coordinates and momenta. The time 
dependent electron-laser pulse interaction is given by 
Hel{t) = Y^kgoca' ^aa'{t)cl^Ck+qa' + h.C. with Vaa'{t) = 
^ / dri(j)oc{riY[A(t) As the coherent phonons 

are generated at g = 0 and the momentum of the light is 
much smaller than the electron momentum, we can set 
g = 0 in all above Hamiltonians. 

We consider parameters in a typical femtosecond 
pump-probe experiment. The pump pulse (central fre¬ 
quency uJq ^ 37bTHz) has a width of r ^ 80/s and 
a relatively broad spectral width of the order of Au ^ 
5 — lOTHz. Such a spectral width is just enough to ex¬ 
cite the lowest energy optical phonon mode whose energy 
is around O ^ bTHz. To ensure that the phonon oscil¬ 
lations are properly resolved in time, the width of the 
pump laser pulse satisfies the condition r « 

The average force driving the coherent phonon oscil¬ 
lations is given by F{t) = —d{He-p){t)/dQ^. Here, 
(...) denotes an ensemble average over eigen states of 
H — He-i{t) time evolving in He-i{t) perturbatively. 
In lines with the authors in ref [38], we assume that 
the electric field is spatially homogeneous and a gaus- 
sian centered around uJq. Thus the electric field product 
E{uj)E{ujF^) is strongly peaked at cjo —f^/2. This leads 
to an expression for the driving force [371138] 

^ -C ^ ( ilnVnm{k)Vml{k) 

iuJo - ^ I ^ - W) 

^ ^ ^ kmn ^ 

^ ^nlYrin{k^^I m W_ 

{^nl + ^ — OJq ig) 



FIG. 1. A cartoon plot showing the toy band structure used to 
illustrate the scattering between superconducting bands close 
to the Fermi level. A quasiparticle is light scattered (solid 
wavy line) from an occupied band state a to an empty state 
in band state b and light scattered again from band state 
b to another band state c. Finally the quasiparticle makes 
a transition back to the band state a by scattering with a 
phonon (dashed-dotted line). The energy scale on the vertical 
axis is of the order of the superconducting gap. 



Here n,m, / are band states, C is an unimportant 
constant, g contributes to the optical absorption, 
and k is the crystal momentum. We have defined 
ujni = 0 Jni{k) = ujn — 0 ^ 1 , where Un is the energy 
of band n with momentum k. The tilde sign above 
the matrix elements denotes the respective quantities 
written in the band basis. In the expression for F{ft), 
we have assumed that the laser frequency is the largest 
energy scale in the problem. Therefore, we have 
chosen to keep the most resonant terms by ignoring 
a third term which has a denominator proportional to 

a) Two band case - analytic result: Our goal in 
this section is to study the phase of the coherent phonon 
oscillations (Fp/^) across Te for a generic two orbital 
model. Our model comprises intraorbital hoppings 
€i{k),€ 2 {k), and the interorbital hopping €i 2 {k) = ruk. 
For analytic simplicity, we choose the two orbitals to 
have the same complex gap order parameter 
This condition will be relaxed in the next section where 
we apply numerics. For the electron phonon couplings, 
we only keep non-zero matrix elements for the coupling 
between the two different orbitals (^') and coupling 
between superconducting particle-hole bands 
As the electronic response to the laser pulse is a very 
fast process, we can assume that the phonon is not 
activated during the laser pulse excitation. Such an 
approximation is easy to justify considering that the 
fast moving electrons have a larger effect on the slow 
moving nuclei than the other way around— an analogue 
of the Born-Oppenheimer approximation in atomic 
physics. In this case, the lowest order effect of the 
electron-phonon coupling is in the form of the driving 
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force, F{t) = -d{He-p){t)/dQ^. 

We now proceed with our calculation of the aver¬ 
age driving force F{t). To bring out the physics essential 
for our discussion, we consider a scattering process 
illustrated in the cartoon in Fig A quasiparticle in 
the state a is scattered by a photon to the empty state 
b above the Fermi level, and then scattered again into 
another empty state c by a second photon. Finally, the 
quasiparticle is scattered back to its original state a 
through a phonon or a series of phonons. We explicitely 
evaluate the matrix element product ^acycb{k)Vi)a{k) for 
such a process so that other similar scattering processes 
can be determined analogously. To do this, we first have 
to perform an unitary transform into the orbital basis 
and then use the formulas described in [39] for tight 
binding matrix elements. We can write the above matrix 
element product as 

LcVeb{k)Vba{k) = - A2) (7) 

x^x_ 

X - A^) +C'sm26»fc(2A,^_)] , 

where f{9k) = -{dmk)^cos'^2dk, x± = ^A^ + S’l, 
(o± = e±{k) + E±{k)^ with the band angle tan20k = 
2mkl{eik — e 2 /c), ^±{k) the band energies, and E±{k) = 

+ €±{k)‘^. From the above expression for the ma¬ 
trix element product, we can separate the most dominant 
contributions from different regions of the Brillouin zone. 
We consider three different cases: (1) contributions from 
momentum space points far away from the fermi sur¬ 
face where e±{k) » |A| > 0, (2) on the Fermi surface 
e+(^) =0 < |A| << e-{k) and finally, (3) on the Fermi 
surface e_(^) = 0 < |A| << 6+(/c). We find that 

iacVcb{k)Vba{k) = ( 8 ) 

+ e±»|A|>0 

fi^k) X I I ^1 + S 2 e) e_ » |A| > 0 = e+ 

e+ >> |A| > 0 = e_, 

where we have defined the effective band energy e = 
and S 2 d = sin20k. Similar expressions can be 
obtained for the other scattering processes. The energy 
denominators appearing in the expression for the driv¬ 
ing force in eq[^ depend quadratically on the energy gap. 
From this, along with the expression for the matrix ele¬ 
ment product (written in eqj^, we arrive at the central 
result of this section — the coherent phonon phase en¬ 
codes the behavior of the superconducting order parame¬ 
ter. For small A, the phase can be written very generally 
as Vph = OLi 0 ^ 2 A(T), where ai and 0^2 are constants 

independent of temperature. As a result, the phase dif¬ 
ference between the superconducting and normal state 


is proportional to the pairing gap. We also addition¬ 
ally conclude that the contribution to the average driv¬ 
ing force from the momentum points far away from the 
Fermi surface is of 0(A/e) smaller than the contribution 
from those close to the Fermi surfaces. However, all the 
regions in the Brillouin zone contribute to the phase of 
the oscillation to the same order. This naturally implies 
that for a significant driving force to be generated, we 
would require the frequency of the phonon mode excited 
5 — lt)THz) to be of the order of the superconducting 
gap. This is a condition that is hard to attain in clas¬ 
sic BCS superconductors, but is comfortably satisfied by 
high Tc Cuprates and iron based superconductors. 

b) Three band case: To further test the above results, 
we consider a more realistic band model that describes 
iron based superconductors and study the pairing sym¬ 
metry dependence. We also examine any signatures that 
can capture the inter-orbital sign change contained in the 
phase of coherent phonon oscillation. To illustrate our 
numerical results and maintain analytical tractability, 
we choose the three band model proposed by Daghofer 
et.al [40|. Fig shows our result for the temperature 
dependence plot of the phase difference F^ — Fat between 
the superconducting and normal states across Tc. The 
phase is a constant above Tc and varies below it due to 
the development of a superconducting gap on the Fermi 
surfaces. For a simple constant s— wave pairing (Fig 

FIG. 2. Plot showing the variation of the phase(rp/i) differ¬ 
ence between superconducting(S) and nornial(N) state as a 
function of temperature across Tc. (Left) Phase as a func¬ 
tion of magnitude of a constant s—wave gap on all the three 
bands. (Right) Phase for different pairing forms of the gap, 
all the same on the three bands. The values of the electron 
phonon coupling is chosen as = 0.4eF for interobital, ^^4 
for xzjyz and ^' /2 for xy intraorbital coupling and the damp¬ 
ing coefficient is chosen as ^ = 0.3eF. The laser and phonon 
frequencies are fixed at 2eF and 0.2eF respectively. 



(Left)), the variation of the phase in the SC state is 
maximum at T = 0 (for small values of the gap) and 
follows a linear dependence on A, as was analytically 
derived in the previous section. However, on increasing 
the magnitude of A, the change in phase develops a 
maximum at a temperature 0 < T < Tc and then falls off 
at T = 0 due to higher order contributions of A. Fig. 
(Right) shows the plot of the phase of the oscillation as a 
function of temperature for different pairing symmetries. 
For the s— wave cases, there is a substantial change in 
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FIG. 3. Plot showing the variation of the real part of the 
driving force iF as a function of temperature (vertical axis) 
and the damping constant g (horizontal axis) in the three 
orbital model of ref |40]. Top row Left (Right): case where the 
signs of the gap on the xz^ yz orbitals is the same (opposite) 
as that on the xy orbital. The color scale represents the real 
part of the driving force iF). Center row (left and right): Cuts 
along different chosen values of g for the corresponding color 
plots above them. Bottom row: The corresponding phases as 
a function of temperature for the p = 0.2 case. The values of 
the electron phonon couplings, laser and phonon frequencies 
are chosen same as in fig|^ 



0.90 0.92 0.94 0.96 0.98 1.00 1.02 1.04 0.90 0.92 0.94 0.96 0.98 1.00 1.02 1.04 

T/r, TIT, 


the phase between T = 0 and T = Tc, whereas for the 
d— wave cases there is little phase change between T = 0 
and T = Tc. In the d— wave scenario, the phase sharply 
plunges on entering into the superconducting state. 

Fig i (top row) shows a color plot of the real part 
of iF as a function of temperature and the absorption 
coefficient g. The left column corresponds to the case 
where the sign of the gap on all the three orbitals is 
the same (+++ case), while that on the right has a 
gap on the xy orbital opposite in sign to that of the xz 
and yz orbitals (++- case). Fig (center row) shows 
cuts corresponding to different values of g for both 
these cases. Clearly, below Tc^ the slope of the real part 
of iF has an opposite sign for the (+++) and (++-) 
cases. More specific to the three orbital model, the real 
part of iF goes through a zero for the (++-) case and. 


therefore, has a tt discontinuity in the phase. On the 
other hand, in the (+++) scenario, the real part of 
iF does not change sign and results in a smooth vari¬ 
ation of phase with temperature (see fig [^(bottom row)). 

To get the physics governing the numerics above, 
we consider the three band model with a definite 
sign of the gap on the xz and yz orbitals (denoted 
by Ai = A and A 2 = A) and an arbitrary gap A 3 
on the xy orbital. We find that for small values of 
A 3 , the driving force on the phonons can be written 
as FiT) = Efc (ai(q + s 5 n(AA 3 ) 4 i(fc)|A 3 (r)|) + 

z (a2(A:) + s5n(AA3)^2(A:)|A3(T)|). Here, di and pi 
are band structure dependent functions which can be 
determined for a given model. For the above model, 
we find that and are both negative. 

This implies that when all the three orbitals have the 
same sign of the gap, the real part of iF(T) is negative. 
On the other hand, if the sign change exists among 
the third orbitals, the denominator becomes zero for a 
critical temperature and results in an observable tt phase 
jump. 

The above results can be applied to investigate the 
pairing symmetries in multi-orbital superconductors. 
Here we specifically discuss its applications to iron-based 
superconductors. Different pairing symmetries, includ¬ 
ing s-wave [4TH46] and d-wave pairing symmetries [47108], 
were proposed for different families of iron-based super¬ 
conductors. Even within the s-wave pairing symmetry, 
there are a variety of possibilities on the internal sign 
changes, including the sign changes between different 
pockets (so called g^ [4Tti44] ) and between bands featured 
by different orbitals (so called orbital-dependent or 
antiphase-g^ [49II52] ). Our results suggest that the phase 
of coherent phonons should have distinct behaviors with 
respect to the antiphase-s^ and d-wave states. In 
particular, as shown in fig[^ if a phase jump can be ob¬ 
served below Tc, it should be a smoking-gun signature 
for the antiphase-5^ state. 

Conclusions'. We have shown that coherent phonon 
oscillations can be an experimental probe of the super¬ 
conducting state and its pairing symmetry. The phase 
of the coherent phonon carries intrinsic information of 
superconducting order parameters and can be applied 
to determine the pairing symmetries in complex multi¬ 
orbital superconductors. 
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